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Inspired by the recent experimental evidence of antiferromagnetism and superconductivity coex- 
istence in heavy fermion CeRhlns, we propose a fully microscopic approach based on the idea of 
real space pairing within the Anderson-Kondo lattice model. We present an overall phase diagram 
incorporating the emergence of a quantum critical point, where Kondo insulating (KI), antiferro- 
magnetic (AF) and superconducting (SC) phases meet. We also obtain the Kondo insulating state 
with totally compensated magnetic moments as the parental state for the emerging SC phases. Fur- 
thermore, the coexistent (AF+SC) phase may contain also a non-trivial spin-triplet gap component 
within the essentially spin-singlet pairing mechanism. 



PACS numbers: 71.27.+a, 74.70.Tx, , 74.20.-z 

The origin of unconventional superconductivity (super- 
fluidity) in strongly correlated electronic [TJ |2] and op- 
tical [21 0] lattice systems is one of the most important 
problems in condensed matter physics, as it concerns go- 
ing beyond one of the most successful theories of XX-th 
century physics - the Bardeen-Cooper-Schrieffer (BCS) 
theory [5] of superconductivity. One of the novel pairing 
mechanisms is based on the Dirac's universal idea of ex- 
change interaction [6 . This idea has a direct application 
in the case of high temperature superconductivity, where 
it takes the form of the t — J model with real space pairing 
[THS]. Here we extend these ideas to the heavy-fermion 
systems [TU] and show that it leads to a unique connection 
between antiferromagnetism (AF) and superconductivity 
(SC) , in accordance with the recent experimental results 
[TT1 [r2] . Also, the evolution from the Kondo insulating 
state (PKI) via pure d-wave SC to the coexistent AF+SC 
phase, leads to highly non-trivial structure of the order 
parameters. 

In this Letter we present a fairly complete phase di- 
agram and thus demonstrate in a straightforward man- 
ner the applicability of the concept of real-space pairing 
to this very important class of strongly correlated quan- 
tum materials. We also show that the system properties 
can evolve from either PKI or antiferromagnetic state to 
the unconventional SC paired states. In this manner, 
we extend the universal meaning of the concept of real- 
space pairing in the strongly correlated systems [101 1131 - 
rT5] , Also, in the PKI state the magnetic moments of 
/ electrons and conduction (c) electrons totally compen- 
sate each other. This last result provides a viable and 
intuitively appealing definition of the Kondo insulators 
and singles out their difference with cither Mott insula- 



tors or intrinsic semiconductors. Furthermore, the paired 
coexistent antiferromagnetic-superconducting (AF+SC) 
state contains also a spin-triplet component within the 
essentially spin-singlet pairing mechanism |16[ I17j . All of 
these results, have been obtained for the first time and 
within a single formal approach. 

The results are discussed starting from introduced by 
us Anderson-Kondo lattice Hamiltonian 10, 18 , in which 
both the Kondo interaction and the residual hybridiza- 
tion processes appear together and are accounted for on 
an equal footing. Subsequently, a modified Gutzwiller 
approach is formulated |19| and combined with the sta- 
tistically consistent renormalized mean-field theory (SC- 
RMFT) [501 IH]- In effect we show that the hybrid (f-c) 
type of real space pairing induced by the Kondo-type in- 
teraction is crucial for description of heavy-fermion SC 
properties even though the /-/ interaction may play the 
dominant role in magnetism of /-electron system in the 
localization limit. We also discuss at the end the role of 
/-/ intersite interaction to demonstrate the stability of 
our solution with respect to this higher-order interaction. 

In the case of cerium compounds, one usually assumes 
that: (i) only the T 7 doublet (4/ 1 ) state of Ce+ 4 -™/ va- 
lency, with the /-level occupancy nj — ¥ 1 in the low- 
temperature dynamics; (ii) this state is strongly hy- 
bridized with the conduction-band (c) states; and (iii) the 
hybridization has either intraatomic or interatomic form. 
We discuss here only the interatomic it leads to 

the stable pairing gap of d-wave character observed al- 
most universally in strongly correlated systems [2]. For- 
mally, in the limit of large intraatomic /-/ Coulomb in- 
teraction of magnitude U the effective Hamiltonian in the 
real-space representation has the following form: 



2 



H — ^ ] \t mn c\ n<J C nc , ^ ] U + 6 ^ i<TC ™°' Cno " ) ^ ; e /^'tcr + ^ ] 



V* V 



U + e f 

£a - «&)M»£<w + h.c) +E J ™ (s< ■ §™ - ^) +E J # • s "'" 

ima im ' 



3 A 



with Jf m ee 2|V im | 2 /(c/ + *f) and j£ ee 
J2mn\ V imVjn\ 2 /(U + e/) 3 . The consecutive terms 
represent the following dynamical processes: the first 
comprises a direct (c-c) hoping in the conduction band, 
as well as the hoping via intermediate /-state. The 
second and the third express, respectively, the bare 
/-electron energy and the spin-flip term. The fourth 
term contains the residual hybridization. The fifth and 
the sixth term represent, respectively, the Kondo (f-c) 
interaction and the Heisenberg (/-/) term (in both 
cases full Dirac-exchange operators are taken). The 
last one represents the Zeeman term for both / and c 
electrons. Here we have projected out completely the 
double occupancies of the / states, what is equivalent 
to assuming that / electrons are strongly correlated. In 
effect the /-electron number n- is replaced by their 
projected counterpart: v ia ee n{ a (l - n{ s ) ee f ia j la 
and Vi ee ^2 a Via- The /-spin operator is defined by 
S, ee (Sf , Sp ee [~f}J lSl l/2& t - %)]. The correspond- 
ing (unprojected) quantities for c-states are n'^ na , 
and s TO . To discuss the system properties a number 
of quantities are first to be calculated self-consistently: 
the /-level and c-band occupancies: nj yC = Ylo-(Ma)> 
the /- and c-magnetic moments: nrif c ee cr(n{^. c ); 
the hybridization correlation and the hybrid pairing 
amplitudes: j a = (f- a c mr7 ) and A CT ee (f ia c mS ), the 
/-/ and c-c hopping correlations x<y = (ficrfja) an d 
£ CT ee (c\ l(J c ma ), as well as the chemical potential /i. 

To calculate the averages of operators (e.g. (fi a c ma )) 
we utilize the extended Gutzwiller scheme [19] . We also 
include the magnitude of the intersite spin-singlet /-/ 
pairing Aff ee (fvffn) and discuss it separately be- 
low. In addition to those, we have to introduce and 
calculate the molecular fields appearing through statis- 
tical consistency conditions so that the self-consistent 
equations for the above quantities coincide with those 
calculated variationally [TS]. Our microscopic model 
is characterized then by the total number of electrons 
n e = (n{ a ) + {n c irJ )\ the magnitude V < of the k- 
dependent hybridization Vk = 8(|V|/W)£k ; bandwidth 
W, and the bare conduction-band dispersion relation 
£k = — (W / /4)(cos(fc x ) + cos(fc y )), (i.e., for exemplary 
two-dimensional lattice, applicable to systems such as 
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FIG. 1: (color online). a): temperature dependence of 
the zero-field magnetic susceptibility for the nonmagnetic 
(parent) Kondo insulator (PKI). The insets show PKI evo- 
lution from the Fermi-liquid regime (FL), to the localized 
moment regime (LM). The upper temperature scale is pro- 
vided for bare c bandwidth W = W 3 K; b): applied field 
h = gusH/ (2W) dependence of magnetic moments and of 
the /-level occupancy (inset). 



CeMIn 5 , with M=Co, Rh, Ir); the position e/ of the 
bare / level; and the magnitude U of the /-/ interac- 
tion. The parameters ej, V, U, and n e are independent; 
also, all the energies are taken in units of W, wherever 
not specified explicitly. 

We discuss first the reference Kondo-compensated 
(PKI) state. In Fig. [T^, we display the temperature 
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dependence of the spin susceptibility \m- The crucial 
point is that the spin susceptibility Xm(T = 0) = 0, as 
one would expect for the totally compensated-moment 
KI state with a small gap or pseudogap. However, in 
distinction to the semiconductors, the system exhibits 
T 2 behavior at low temperature charactering the Fermi- 
liquid state. Also, the x m (T) exhibits a maximum for 
T max = 50+250.K (when W is in the range - 0.1-0.5eV) 
followed by the Curie- Weiss behavior with a large (and of 
AF sign) paramagnetic Curie temperature in the range 
Q = 100 — 150K. These features, together with the uni- 
versal scaling law XmP = const., with p being the sys- 
tem resistivity, proposed and tested earlier [22], provide 
us with the confidence that our method of approach de- 
scribes properly the universal trend of the experimental 
data for KI systems [35]. Furthermore, we have shown 
Fig. [TJd the applied magnetic- field dependence of the /, c, 
and the total magnetic moments at T = 0. The total mo- 
ment remains completely compensated in low fields and 
this is a true Kondo-lattice state, broken at the meta- 
magnetic point fXBH m /W ~ 0.02. The magnetization 
curve above metamagnetic critical field H m shows ap- 
proximately linear behavior, m t ~ H, as marked. In the 
inset we show field dependence of the /-level occupancy 
(valency). Note that even though the ground state is in- 
sulating, the valency deviation from Ce 3+ configuration 
is 1 — rif = 8 ~ 0.2 i.e., is not integer (n/ = 1 case 
corresponds to the Ce 3+ localized-moment state and is 
reached only as |V| — > 0). 

We turn now to the analysis of magnetic and super- 
conducting phases appearing in the metallic state, i.e., 
with the partial filling of orbitals, n e < 2. In Fig. [2] we 
draw an overall phase diagram on the plane occupancy 
n e - hybridization magnitude \V\. This diagram con- 
tains two coexisting AF+SC phases in the heavy-fermi- 
liquid regime (n e -+ 2), as well as AF and SC phases. 
In the small-hybridization limit we encounter two fer- 
romagnetic phases (WFM and SFM). In insets of the 
magnetic part we characterize the phases by the corre- 
sponding shapes of the quasiparticle densities of states 
and the position of the chemical potential. Also, in the 
lowest right inset we show the evolution of the correlated 
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and (Aj.) c = A it = (/i[.c mt ) c . They are of equal am- 
plitude and of opposite sign only when the supercon- 
ducting state is of pure spin-singlet nature. The phase 
diagram drawn here along the line n e ~ 1.95 reflects 
the sequence of phases observed recently in CeRhIn 5 
(AF— !>AF+SC^SC) with the increasing pressure [TT]. In 
general, as in the case of high temperature superconduc- 
tors, pure SC state evolves with the diminishing electron 
count from the insulating, here nonmagnetic PKI state. 
Also, in distinction to the high-T c systems, we have here 
two AF kinetic-exchange interactions: the Kondo f-c and 
/-/ type, both arising from the virtual inter-band c f 
hopping. To characterize the phases marked in Fig. [2] 
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FIG. 2: (color online). Overall phase diagram on the plane 
number of electrons n e = (n{) + (nf) - interatomic hybridiza- 
tion magnitude |V|, comprising strong ferromagnetic (SFM), 
weak ferromagnetic (WFM), coexistent (AF+SC), antiferro- 
magnetic (AF), superconducting (SC), and the parent Kondo 
insulating (PKI) states. Inset: evolution of the SC state from 
PKI; the spin components of (A) c the spin-singlet gap param- 
eter: (Af) c = A n = (/itc m 4.) c and (A4,} c = Aj.-f = {fnc m f) c 
are shown. The solid dot marks the AF metal - Kondo in- 
sulator (PKI) quantum critical point (QCP). The parameters 
are U = 2, e/ = —0.85. The representative density-of-states 
curves define the normal-phase distinctions. 
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FIG. 3: (color online). Magnetic moments as a function of 
|V|. The left panel presents the uniform component m FM 
of magnetization in the specified phases. The top right panel 
shows the /-level occupancy n/ as a function of the hybridiza- 
tion and the bottom represents the staggered component of 
magnetization, m AF . The magnitude of the uniform compo- 
nents of magnetic moments was multiplied by the factor 10 
in the AF+SC and SC phases. 



we have drawn in Fig. [3] evolution of the calculated 
uniform (m FM ) and the staggered (m AF ) components 
of the magnetic moment, defined through the relation 
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»V(c)<<t = W<T')o = |("/( c ) + -'"/(c) 
with Q = (7r, 7t) being the AF-ordering wave vector and 
Ri - the lattice-site position. Fig. [3] conveys our next 
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FIG. 4: (color online), a) illustration of the two different AF 
orders as a consequence of both the /-/ (oc J H ) and Kondo- 
type f-c (oc J K ) spin-spin interactions; b) orientation of the 
/-moments (blue) and the c-moments (red) in FM and AF 
phases, respectively; c) different components of the hybrid (/- 
c) superconducting gap in the AF+SC phases, which compose 
the spin-singlet Af = (A^ + Af )/2 and the spin-triplet A^ = 
(A^ — Act)/2 components of the gap. 
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FIG. 5: (color online). Temperature dependence of the 
correlated-gap components, as well as of m AF and m FM (in- 
sets). The parameters: n e = 1.95, \V\ = 0.28. The arrows in 
the AF+SC phases label Af' T components with a =\ or J,. 

principal messages coming out from our fully microscopic 
theory. First, the uniform moment component is practi- 
cally compensated in the AF+SC phase (the exemplary 
numerical values are listed in table 1 of Supplementary 
Material). Second, there appears a very small staggered 
moment in the second to-the-right AF+SC state. This 
may explain, why it may be difficult to measure it in 
the coexistent phase as discussed by Knebel et al. [TT] . 
Third, most of the transition lines mark discontinuities 
and thus, their first-order nature as a function of | | 
(emulating the changing pressure applied to the system) . 



Fourth, the coexistent AF+SC phases border with AF 
and SC states and have a very small uniform magnetic 
moment. In general, the diagonal (spin-spin) correlations 
lead to the Kondo compensation and the off-diagonal 
(pairing) correlation lead to the emergence of supercon- 
ducting state. Finally, the valence of Ce +4_n/ evolves 
systematically from that for heavy- fermion rif = 1 — 6, 
with 5 <C 1, to that characterizing the fluctuation- valence 
regime (rif < 0.75), with the increasing hybridization. 

Before characterizing the temperature evolution of the 
gap components, we sketch in Fig. [4] both configura- 
tions of the component moments in AF and FM phases, 
as well as define different gap components (c), that we 
have to introduce in order to describe the coexistent 
AF+SC phases. Such multiplicity of the gaps is en- 
riched when we include also the d-wave gap coming from 
the /-/ pairing. It is important to note that the spin- 
triplet component appears only in the coexistent AF+SC 
phase |17j . as there are majority and minority spins 
then, e.g. on sublattice A we have that (h^) > (n^) 
and similarly for the c electrons; hence A^ ^ A^ . In 
Fig. [5] we show that all the gaps define still a sin- 
gle transition temperature for given phase. One should 
note the sequence AF+SC— >SC— >-AF— ^paramagnet. 
AF— >paramagnet transition temperature is of an order of 
magnitude higher than that for AF+SC— >SC transition. 
The /-/ pairing coming from the fourth-order contribu- 
tion of magnitude jffj\ = 12V 4 /(e/ + U) 3 is also taken 
into account in this panel and, in distinction to high- 
T c systems, here the hybrid (f-c) pairing contribution is 
dominant. In all the phases we have that rif ~ 0.9 + 0.8. 

To recapitulate, the strong correlations, crucial for the 
appearance of both the AF and the Kondo-compensated 
states, are responsible to an equal extent for the emer- 
gence of both the coexistent AF+SC and the pure SC 
phases. The sequence of the phases AF— >■ AF+SC— >SC 
appears with the increasing V/tf ratio, as observed in 
quasi-two-dimensional Cefthln5 [TTJ [T21 [23] • Our real- 
space pairing thus accommodates in a natural manner 
into the spin correlations in AF phase (cf. Fig. [4^). The 
interesting feature is also the presence of a very weak 
staggered component of magnetic moment in the second 
AF+SC phase (cf. Fig. [3]). The principal result is also 
the strong mutual /- and c-moment compensation, either 
complete or partial, depending on the phase. Further- 
more, the phenomenon assisting the uniform moment is 
the appearance of the spin-triplet gap component in the 
coexistent AF+SC phases. 
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